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Abstract 



This paper considers random cellular systems in one, two, or three dimensions. We study the carrier- 
to-interference ratio (CIR) and the carrier-to-interference-plus-noise ratio (CINR) at the mobile-station 

G\ ■ 

(MS), by using the so-called shotgun cellular system (SCS) model for the cellular system. In a SCS, 
the base-stations (BS) are placed randomly according to a non-homogeneous Poisson point process. For 



the SCS, analytical expressions for the tail probabilities of CIR and CINR are obtained. These results 
are generalized for an SCS affected by random shadow fading with any general distribution, for any 
arbitrary path-loss models, and for the MS located anywhere in the SCS. Further, for a special case of 
the SCS, namely, the homogeneous SCS, it is shown that the shadow fading has no effect on the CIR, 
and the effect of random transmission gains, power control, multiple channel reuse groups are analyzed. 
The techniques developed for the analysis of the CIR and the CINR holds applications beyond cellular 
networks and can be used to study several important scenarios in both uplink and downlink for cognitive 
radios, femtocell networks and other heterogeneous and multi-tier networks. 

Index Terms 



Cellular Radio, Co-channel Interference, Fading channels, Stochastic ordering. 
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I. Introduction 

The modern cellular communication network is a complex overlay of heterogeneous networks 
such as macrocells, microcells, picocells, and femtocells. The base station (BS) deployment 
for these network can be planned, unplanned, or uncoordinated. Even when planned, the base 
station (BS) placement in a region deviates from a regular hexagonal grid due to site- acquisition 
difficulties, variable traffic load, and terrain. The coexistence of heterogeneous networks have 
further added to these deviations. As a result, the BS distribution appears increasingly irregular 
as the BS density grows and is outside standard performance analysis. 

Two approaches of modeling have been widely adopted in the literature. At one end, the BSs 
are located at the centers of regular hexagonal cells to form an ideal hexagonal cellular system. 
At the other end, the BS deployments are modeled according to a Poisson point process which 
we denote a shotgun cellular system (SCS). In (TJ, the author makes a connection between 
these two models on a homogenous two dimensional (2-D) plane. It is shown that the carrier-to- 
interference ratio, ( j), of the SCS lower bounds that of the ideal hexagonal cellular system and 
moreover, the two models converge in the strong shadow fading regime. Since the BS deployment 
in the practical celluar system lies somewhere in between these two extremes, it is important 
to completely understand the performance of the SCS. Such an analysis holds significance in 
areas beyond cellular networks. Wireless LANs, cognitive radios, and ad-hoc networks are also 
characterized by irregular deployment of the BSs [|2ll- [fT3l . It is emphasized that, although the 
deployment of BSs in practice is not random, such a study is useful because it allows a thorough 
theoretical understanding of many important effects in the strong shadowing regime. 

Recent work and our contributions: A Poisson point process has been adopted in the literature 
for the locations of nodes in the study of sensor networks, ad hoc networks and other uncoor- 
dinated and decentralized networks. In the case of ad-hoc networks, bounds on the transmission 
capacity are derived in several different contexts Q, flU, [fTTTl . ffT2ll . lfl4l - lfT8i Finding the 
optimal bandwidth partitioning in uncoordinated wireless networks is considered in |[T9l . Similar 
outage probability analysis in ad-hoc packet radio networks is considered in COfl , [1211 . 

An underlying assumption in all the previous work is that the density of transmitters is 
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constant throughout the cellular region, i.e. the Poisson point process is homogeneous. Such 
a model does not appropriately represent practical cellular networks where the BS distribution 
is irregular. In this paper, this scenario is incorporated by modeling the distribution of BSs by a 
non-homogeneous Poisson point process. Moreover, the region of interest need not be restricted 
to IR 2 as in prior work, and may be IR 1 or IR 3 . Furthermore, the performance dependence on the 
MS location within the non-homogeneous cellular region is also characterized. Note that only 
the instantaneous j at the MS are captured here. Handoff features and other dynamics are out 
of the scope of this work. 

Finally, most research restricts the interference analysis to popular fading models like log- 
normal, Rayleigh, and Rician distributions and a propagation model that follows the inverse 
power law path-loss. Here, results are obtained that hold for any general fading distribution and 
arbitrary propagation models. This facilitates more accurately modeling the wireless network for 
analysing the signal quality at the desired receiver. 

The following are the main contributions of this paper. 

• A semi-analytical expression for the j and distribution at the MS in a SCS with the 
propagation loss modeled by an inverse power law. (Sections IIII-AI and UVT) 

• An accurate closed form approximation to the SCS j distribution. (Section IIII-B|) 

• The SCS y and distributions under shadow fading with any probability distribution 
(Section [V]), and for any arbitrary propagation loss model (Section IVIl) . 

Finally, in Section rvTTl applications of the results listed above in specific wireless communication 
networks are briefly described. In the following section, the system model is described and the 
performance metric of interest is defined. Further, the model is reduced to an equivalent Poisson 
point process in 1-D that is the basis of analysis in subsequent sections. 

II. System Model 

This section describes the various elements used to model the cellular system, namely, the 
BS layout, the radio environment, and the performance metrics of interest. In the following 
subsection, we define SCS and introduce the terms associated with the BS layout in a SCS. 
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A. BS Layout 

Definition 1. The Shotgun Cellular System (SCS) is a model for the cellular system in which 
the BSs are placed in a given /-dimensional plane (typically I — 1,2, and 3) according to a non- 
homogeneous Poisson point process on M. 1 . The intensity function of the Poisson point process 
is called the BS density function in the context of the SCS. 

In a practical cellular system, for places where the mobile user population is very small, the 
BSs are typically deployed along the highways. Such a cellular system can be modeled by a one- 
dimensional (1-D) SCS. . In the 1-D SCS, the BSs are placed according to a non-homogeneous 
Poisson point process with a BS density function A(r) > 0, where r > is the distance of the 
BS from the origin on R 1 , such that the probability that a BS is placed in an interval (r, r + h) of 
infinitesimal length h, is X(r)h. The formal definition of the 1-D SCS incorporating the definition 
of the non-homogeneous Poisson point process is provided below. 

Definition 2. A 1-D SCS with an average BS density function A(r), where r > is the separation 
of the BS from the origin, is a counting process, {iV (r) , r > 0}, where N (r) is the number 
of BSs upto a distance r from the origin, which satisfies 

N(0) = 0, 

{N (r) , r > 0} has independent increments, 
P({1 BS in (r,r + h)}) = X(r)h + o(h), and 
P ({0 BS in (r, r + h)}) = 1 - A (r) h + o (h). [|2l 

Similarly, a two-dimensional SCS (2-D SCS) is a model for the planar deployment of BSs. 
Let every point on this 2-D plane be described by the polar coordinates (r, 9). In the 2-D SCS, 
the BSs are placed according to a non-homogeneous 2-D Poisson point process with a density 
function A(r, V r > 0, < 6> < 2tc, which is the density of BSs about the origin in M 2 at 
the point (r, 9), so that the probability that there exists a BS in an infinitisimal area dA = rdrdO 
about the point (r, 9) is A (r, 9) dA. 

Further, a three-dimensional SCS (3-D SCS) models the BS deployments in metropolitan areas, 
or wireless LANs (WLAN) in residential areas. If every point on the 3-D plane is described by 
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the spherical coordinates (r, 9, 0), then the BS placement in the 3-D SCS takes place according to 
a non-homogeneous 3-D Poisson point process with a density function A (r, 9, 0) , V r > 0, < 
9 < n, < < 2ir, which is the density of BSs about the origin in 3-D at the point (r, 9, 0), 
such that the probability that there exists a BS in an infinitisimal volume dV = r 2 sin9drd9d<p 
about the point (r, 9, 0) is A (r, 9, 0) dV. 

A 1-D SCS in which A (r) is a constant everywhere on M 1 is termed as a homogeneous 
1-D SCS. Similarly, a 2-D/3-D SCS in which the respective density functions are a constant 
everywhere are termed as homogeneous 2-D/3-D SCS. 

B. Radio Environment 

In the SCS, all BSs are assumed to have identical transmission powers and antenna gains, and 
omni-directional antennas. The signal from the BS propagates through the wireless channel where 
it undergoes a deterministic path-loss, a random shadow fading and is affected by background 
noise. The path-loss is modeled by an inverse power law, and is proportional to R~ £ , where 
e is the path-loss exponent, and R is the distance from the BS. The shadow fading factor is 
often modeled as a random variable with a log-normal distribution [|23l , E4l . The thermal and 
background noise in the system is given by a constant noise power N. Overall, the signal power 
at a distance R from the BS, accounting for both the path-loss and the shadow fading, is 

P = K^R~ £ , (1) 

where K captures the transmission power and the antenna gain of the BS, and ^ is the random 
shadow fading factor. 

C. Performance Metric 

In this paper, we focus on the downlink performance of the SCS. In other words, we are 
concerned with the signal quality at a mobile- station (MS) within the SCS. The MS is without 
loss of generality located at the origin of the plane in which the SCS is defined, unless specified 
otherwise. The MS chooses to communicate with the BS that corresponds to the strongest 
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received signal power. This BS is referred to as the "serving BS", and all the other BSs in 
the SCS are collectively called the "interfering BSs". Consequently, the signal quality at the 
MS is defined as the ratio of the received power from the serving BS (C or Pg) to the sum 
of the total interference power (/ or Pj), and the noise power (N), and is called the carrier- 
to-interference-plus-noise ratio (jT]v). In an interference limited system, I ^> N and the signal 
quality is the carrier-to-interference ratio (j). 

Using (Jlj), the y and at the MS are expressed as 

C Kg^gR- s £ C Kg^gR- s £ 

i TZx K ^i R r' m i + n E^vMT + iv' {) 

where subscript S denotes the serving BS and subscript i indexes the interfering BSs; Kg, {K i } c *L 1 
are independent and identically distributed (i.i.d) random variables; R s , are random 

variables that come from underlying Poisson point process that govern the BS placement; 
^g, {^i}^ are independent and identically distributed (i.i.d) random variables. Hence, j and 
are random variables, and can be characterized by a probability density function (p.d.f.), 
cumulative density function (c.d.f.) or the tail probability. The tail probability of is given 
by P ({j+n > ^})' arK * ls tne probability that a MS in the SCS has a signal quality of at least 
r], 7] > 0. Next, we characterize the tail probability of the y at the MS in a SCS. 

III. j CHARACTERIZATION IN THE SCS 

Here, the transmission power and the antenna gains of all the BSs in the SCS are assumed 
to be constant (say, K). Also, the shadow fading factors are assumed to be a constant. Hence, 
from the expression for j in (j2J), the BS closest to the MS is the serving BS and the expression 
for j can be equivalently written as 

C KR^ £ 



-e 



(3) 



/ TZ2KR 

where R 1 < R 2 < R3 • • • are distances between the BSs and the MS in a non-decreasing order. 

In Section ITlI-Ai we first obtain the tail probability of j for a MS at the origin of the 1-D 
SCS and illustrate the procedure to obtain the same when the MS is not located at the origin 
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(refer Section ITlI-A II) . Then, we obtain similar results for the 2-D and 3-D SCS by transforming 
them to an equivalent 1-D SCS (Section [III-A2I) . In Section Ull-Bl we focus on studying the j 



and 



c 

I+N 



in the homogeneous 1-D, 2-D and 3-D SCS. 



A. Tail probability of j 

1) 1-D SCS: To obtain the tail probability of the j at the MS in the 1-D SCS, we first 
present a theorem that derives the characteristic function of (~) 1 at the MS. In order to state 
the theorem, we need the following definition. 

Definition 3. A 1-D BS density function A(r), r > is said to be stable with respect to the 
path loss exponent e if JJ^ X(s)ds < oo and ^lim f^ ri X(r)r~ £ dr < oo, V n > 0. 

Theorem 1. In the 1-D SCS with a BS density function A(r) that is stable with respect to e, 
(a) the characteristic function of Pi conditioned on Ri is given by 

{iuK) 



Pi\Rx {^\ r h 



v n=l 



111 



/ A (r) r n£ dr J , and 

J r=r 1 J 



(4) 



(b) the characteristic function of (y) is given by 

^o)-i(w) = E Rl 



■PARi 







(- 







exp x j (exp (iuju £ ) — l) A (uRi) du 

where [■] is the expectation operator with respect to the random variable R\. 

Proof: See Appendix |Bj 
The tail probability of j may be directly obtained from the characteristic function and is 



(5) 
(6) 



Hy ' '/ 




,-i [CO 



1— expl 



duj 
2 



T] = 0. 



(7) 



The above has been previously proved (See [251 Eq. (9)]). 
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Now, consider the problem of finding the tail probability of the j at the MS in the 1-D SCS, 
but not located at the origin of the 1-D plane. Notice that A(r) has been defined as a density 
function w.r.t. the distance from the origin. When the BS density is defined for the entire 1-D 
plane (d(x), — oo < x < oo), i.e. for either sides of the origin, the corresponding A(r) is 

A(r) = d(r) + d(-r), r > 0, (8) 

which is the sum of the densities of all the points that are at a distance r from the origin in the 
1-D plane. This corresponds to x = r and x = — r, respectively. The following lemma helps in 
finding the BS density as a function of distance from a MS at any point y on the 1-D plane. 

Lemma 1. Given an arbitrary density function d(x), — oo < x < oo defined for the entire 1-D 
plane, and a MS at y on the 1-D plane, 

A(r) = d(y-r) + d(y + r), r > (9) 

is the density function as a function of the distance from the MS. 

Proof: As in ([8]), A(r) is the sum of the densities of all the points that are at a distance r 
from the MS, and this corresponds to the points y + r and y — r in the 1-D plane. Hence, A(r) 
is given by (j9]). ■ 
Now, given d(x), — oo < x < oo, finding the tail probability of j at the MS located at an 
aribitrary point y on the 1-D plane is a three step procedure. First, use Lemma [T] to find the 
equivalent A(r). Second, use A(r) in Theorem \T\ to obtain $/ C n-i (uj), and thirdly, evaluate this 
function in (jTJ) to obtain the desired result. The three step procedure now allows us to obtain the 
tail probability of the instantaneous j at the MS which is anywhere along the 1-D SCS. Next, 
we show an interesting property of the j at the MS in the 1-D SCS, for which we require the 
following definition. 

Definition 4. Let X and Y be two random variables such that P ({X > x}) < P ({Y > x}) , 

V x G (— oo, oo) , then X is said to be smaller than Y in the usual stochastic order and is given 
by X < st Y. Also, X = st Y means P ({X > x}) = P ({Y > x}) , V x E (-oo, oo). fl26]| 



Corollary 1. The p.d.f. (c.d.f and tail probability) ofj at the MS in the 1-D SCS with BS density 
function A(r) is the same as that in 1-D SCSs with BS density functions -A(-), V a > 0, i.e., 



c\ c I 

l| A (r) °< l\^X(ry 



Proof: Let {-Rfcl^Li correspond to the 1-D SCS with BS density function A(r).Then, since 
the ordered base station locations Rk's, are determined by inter-base station distances, it follows 



that Qr I 



W (aRi)- 



(6) Ji,, 



A(r) St ESLaK) 



, where (a) is obtained by multiplying 

L*(£) 



the numerator and denominator by a~ £ , a > in ([3]) and (6) follows from Appendix lAl-Lemma 
E3 Further, {i4)d in (6) correspond to 1-D SCS with BS density ±A ( r -) , a > 0. U 

As a result, {^A (^) , r > 0, a > 0} forms a parametric family of BS density functions 
such that the 1-D SCSs corresponding to them have the same y at the MS. In other words, 
appropriately scaling the BS density function will not change the p.d.f. of j. Next, we analyze 
the tail probability of j at the MS in the origin of 2-D and 3-D SCS. 

2) 2-D and 3-D SCS: Notice from the expression of y in (j3]), that the distance of the BS 
from the MS is sufficient, and the additional details about the orientation of the BS in the 2-D or 
3-D plane is irrelavant. Hence, for the study of the y, the 2-D and 3-D SCS can be transformed 
to an equivalent 1-D SCS using the following lemma. 

Lemma 2. (a) The distances of the BSs from the origin, r, in the 2-D SCS with a BS density 
function A(r, 9) forms a Poisson point process with density function A(r) = J^ Q A(r, 9)rd9, r > 
0, provided the integral exists. Hence, for the y analysis, the above 2-D SCS is equivalent to 
1-D SCS with BS density function A(r). 

(b) The distances of the BSs from the origin, r, in the 3-D SCS with a BS density function 
A(r, 9, (j)) forms a Poisson process with density function A(r) = J^ =Q J^ Q A(r, 9, (f))r 2 sin9d(j)d9, 
r > 0, provided the integral exists. Hence, for the y analysis, the above 3-D SCS is equivalent 
to 1-D SCS with BS density function A(r). 

Proof: Lemmata) and|2|&) are a corollary of the Mapping theorem for Poisson processes 
in [|27l Page 18], which says that the projections of a higher-dimensional Poisson processes are 
also Poisson processes. ■ 



Once the 2-D and 3-D SCSs have been translated to an equivalent 1-D SCS, the problem of 
finding the tail probability of j is solved using the results in Section IHI-All Now, we focus on 
further studying the characteristics of j at the origin of a homogeneous l-D SCS, Z = 1,2,3. 

B. Homogeneous l-D SCS 

An important features of the homogeneous l-D SCSs is illustrated by the following lemma. 

Corollary 2. A homogeneous l-D SCS with a constant BS density A over the entire space is 
equivalent to the l-D SCS with a BS density function A;(r) = \ bir , V r > 0, where bi = 2, 
bi = 2n, 6 3 = 4-7T. 

Proof: For the homogeneous l-D SCS, the result follows from ([8]) with d(x) = A , V x. 
For the homogeneous 2-D SCS and 3-D SCS, the result is obtained by evaluating the integrals 
in Lemma [2{a) and W(b), respectively. ■ 
Notice in Corollary [2] that the BS density function of the equivalent l-D SCS of the homoge- 
neous l-D SCS is parameterized by Ao and I. Consider setting up Monte-Carlo simulations for 
the homogeneous l-D SCS for comparing the tail probability of j as a function of A . Figure 
Q] is the result of such an experiment and it shows that the tail probability does not vary as a 
function of A . The following corollary warrants this empirical result. 

Corollary 3. In a homogeneous l-D SCS with a constant BS density X , if the path-loss exponent 
satisfies e > I, 

(a) the characteristic function of Pi conditioned on R\ is given by 

$p iIRi (u\n) = exp -j Fi 1 - I iuKrA ^ , and (10) 

(b) the characteristic function of (j) is given by 

1 



®Pi\Ri 



P~s 



Ri 



(11) 



1 F 1 (-L-l-L-iu;)' 

where E Rl is the expectation w.r.t. R lt and iFx (■;■;■) is called the confluent hyper geometric 
function of the first kind (Refer ll28l Chapter 1, Section 1.7] for the definition). 
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Proof: Using Corollary [2l the density function of the equivalent 1-D SCS is obtained, A(r) = 
Xobir 1 " 1 . Next, by evaluating the integral in Theorem [T^a) , ( JlO|) is obtained. Finally, using Fact 
|2]in Appendix|Al compute the p.d.f. of Ri to be /^(n) = Xokr^ 1 exp (-^f 1 ^) , < r\ < oo, 
using which the expectation in ffTT]) can be evaluated to obtain the final result. ■ 
The significance of Corollary [3] is presented in the following two remarks. 

Remark 1. The characteristic function of the ( j) 1 does not depend on A , and hence the tail 
probability of j at a MS in a homogeneous /-D SCS does not depend on A . 

Remark 2. The expression for the characteristic function of ( j) 1 for a homogeneous 2-D and 
3-D SCS is same as of a homogeneous 1-D SCS with path-loss exponents | and |, respectively. 

Remark Q] proves why the curves corresponding to the homogeneous 1-D, 2-D and 3-D SCSs 
in Figure [Dare constant values. Remark [2] helps build an intuition of why the homogeneous 1-D 
SCS has a higher tail probability of j than homogeneous 2-D and 3-D SCSs; Figure Q] now 
corresponds to comparing the tail probabilities of j in a homogeneous 1-D SCS with path-loss 
exponents e, |, and |, respectively. As the path-loss exponent decreases, the BSs farther away 
from the MS have a greater contribution to the total interference power at the MS, and this leads 
to a poorer ~ at the MS and a smaller tail probability. 

In Figure [21 the tail probability of j is plotted against e; the squares (□) and the pluses 
(+) represent the values computed using (j9]) and Monte-Carlo simulations, respectively. The 
advantage of using ()9]) becomes evident when the path-loss exponent is arbitrarily close to 2 for 
the homogeneous 2-D SCS case, where the convergence of the Monte Carlo simulations is very 
slow and takes a large number of trials. An important corollary of Remark [2] is discussed below. 

Corollary 4. For a homogeneous l-D SCS, I = 1, 2, 3 ■ ■ ■, where the path-loss exponent is e, the 
tail probability of j at the MS is given by 

F ({t > ^j) = p ({t >1 }) xr, ~'' Vr? - 1) £>l (12) 

= £|7T«, V?7> 1, £>l, (13) 
where KLl is a constant parameterized by f. 
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Proof: In HI, we have shown that 

¥ {{ < 7 >V }) = F UJ >1 }) xr?_f ' V7? - 1, £>2 ' (14) 

for a homogeneous 2-D SCS. From Remark [2l ( fl4]l holds for all homogeneous 1-D SCS with 
path-loss exponent | and therefore, for all homogeneous 3-D SCSs with path-loss exponent y. 
Hence, (fT2|) hold true, and (fT3j) is obtained by noting that the characteristic function of ( ~) 1 
is a function of | and soP({j>l})isa constant parameterized by f. ■ 
Recall that /C| can be obtained by using ( fTTl) and evaluating the integral in (jTJ) with 77 = 1. 
Note that j is a non-negative random variable with a support of [0, 00), and surprisingly, its tail 
probability has such a simple form as given by (|T3|) in the region [1, 00). Next, we define the 
so-called few BS approximation and derive closed form expressions for the tail probability of — 
at MS in a homogeneous l-D SCS for both the regions [0, 1) and [1, 00). 

Definition 5. The few BS approximation corresponds to modeling the total interference power 
at the MS in a SCS as the sum of the contributions from the strongest few interfering BSs and 
an ensemble average of the contributions of the rest of the interfering BSs. 

Recall from (j3]) that Pj = KR^ e , where {Ri}^ is the set of distances of BSs arranged 
in the ascending order of their separation from the MS. The arrangement also corresponds to 
the descending order of their contribution to Pj, due to path-loss. In the few BS approximation, 
Pj is approximated by P l (k) = £-=2 KR~ £ + E [EZk+i KR T £ \ R k] , for some k, where E [■] 
is the expectation operator and corresponds to the ensemble average of the contributions of BSs 
beyond Rk- The j at the MS obtained by the few BS approximation is denoted by j^-. 

Corollary 5. For a homogeneous l-D SCS, with BS density Ao and e > I, for k=l,2,3,- ■ ■ 

X biKR l r £ 



E 



E KR i 

J=k+1 



Rk 



(15) 



Proof: Firstly, use Corollary [2] to reduce the SCS to the equivalent l-D SCS with BS density 
function A(r) = Ao^/r^ 1 , V r > 0. Next, given k, use Fact[3]in Appendix |A] to divide the [0, 00) 
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to non-overlapping intervals [0, Rk] and (J4, oo), such that the events in them are independent of 
each other. Next, from Definition [2l the probability that there exists a BS in (r, r + dr) is \(r)dr. 
Hence, the ensemble average of the contribution a Pj due to BSs beyond Rk is the weighted 
sum of the contributions to Pj due to a BS in (r, r + dr) for every r e oo) weighted by 
A(r)dr. Hence, £ [E*= fc +i #fc] = K\ Q bir' £+l ~ l dr and flB]) is obtained. ■ 



Now, we present the theorem to obtain the tail probability of jr- = j- 



in closed form. 



k=2 



Theorem 2. In a homogeneous l-D SCS with BS density Xq and path-loss exponent e, satisfying 



e > I, the tail probability of ^ at the MS is given by 



P| <y 2 >v 



7] eC|, 



V > 1 



1 



-«(»») (i + 14(77)) + rj~^Ds (77), 77 < 1 



(16) 



w/zere C« = G(0) and D« (77) = G(m(?7)) wfcfc G(a) = JZ a ve ~ V i 1 + v (j ~ X ) 

« (17) = (f-i) (^-i N 



cfo, and 



Proof: Here we provide an outline of the proof. First, using Corollary [51 



/ 2 Pr(2) 



with 



Pi (2) = Ki? 2 6 (l + 717 #2) • Next > notice that the event j £ > 77} is equivalent to the joint event 



Ri < R 2 , -Ri < 
where 



min R 2 , 



K 



and thus, P 



K 



T 2 >V 



K 
K 



P 



R x < min ( R 2 , ( ^ 



T7 > 1 



lxu(r]) \ ; 



, 77 < 1, # 2 < 



lxu(ij) \ 1 



Finally, (JTHJ) is obtained using the joint p.d.f., f Rl>R2 (n,r 2 ) = (A fy) 2 (r 1 r 2 ) /_1 exp (-^r 2 ) , 
< 7*1 < r 2 < 00, which results from Fact Q] and Fact [2] of Appendix lAl ■ 
Notice that P ({^ > 77}) = ^P ({£ > f° r ?7 > 1, and Figure |2] shows that the few 
BS approximation (•) closely follows the actual behavior (□). Figure |3] shows the comparison 
of the tail probabilities of j (computed using the characteristic function of (y) 1 ) and j- for 
a homogeneous 2-D SCS with path-loss exponent 4. Notice that the gap between the two tail 
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probability curves is negligible in the region r\ G [0, 1], and further, both the curves are straight 
lines parallel to each other in the region r\ E [1, oo), when the tail probability is plotted against 77, 
both in the logarithmic scale. This shows that the few BS approximation, is a good approximation 

for the actual j and also characterizes the signal quality in closed form. Next, we use the results 

c 

I+N 



developed so far, to study the characteristics of at the MS in a SCS 



IV. CHARACTERIZATION IN A SCS 

Without loss of generality, we study the at the MS in the 1-D SCS. We first obtain the 
tail probability of using the characteristic function of (j£v) 1 derived in the following 
corollary. 

Corollary 6. In the 1-D SCS with a BS density function A(r) that is stable with respect 
to e, the characteristic function of the sum of the total interference power (Pi) and noise 
power (N) conditioned on Ri is $p /+ jy|i? 1 (w| ri) = exp (iuN) x & Pl \ Rl (u)\ri), and the 



characteristic function of (7^7) 1 is ®r_c_\- 1 ( w ) = exp (ip^Nj x &Pi\Ri (7^: 
where $p 7 |Hi ( w l r i) 1S given by (j4j). 



Ri 



Proof: The expressions for ^p i+ n\R! (^I r i) an d < £/_^\- 1 ( w ) follow directly from the 



definition of characteristic function, where TV is a constant. 

c 

I+N 



Further, the tail probability of -J^rr is 



' T] = 0. 

and is obtained by substituting — with in (jTJ). Next, an interesting property of the at 
the MS in a homogeneous l-D SCS is given in the following corollary. 

Corollary 7. If the at the MS in a homogeneous l-D SCS is specified by (A , e, K, N) where 
A is the constant BS density, e is the path-loss exponent, K is transmission power and antenna 
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gain of each BS, and N is the constant noise power, then, 



C 



Proof: 



c 



I + N 



(a) 



C 



"St 



(X ,e,K,N) 



(ggiT 



I + N 



l,e,l,NX 7 K- 1 



(18) 



I+N\(\ , £ ,K,N) YZtaicJUr 



N 

' {Ka e ) 



(b) 



( i 



-st 



Mr) 



where a 



a 1 V a / 



Aq , K —1, N = N (Ka e )~ , (a) is obtained by expressing in terms of the equivalent 1-D 
SCS with A; (r) = Aofyr* -1 , r > 0, and multiplying numerator and denominator with (ifa e ) _1 , 
(6) follows from Corollary [Q and finally, (HHJ) is obtained by noting that the 1-D SCS with BS 
density function -A; (^) in (6) corresponds to a homogeneous /-D SCS with BS density 1. ■ 
So, it is sufficient to analyze a homogeneous l-D SCS with BS density Ao = 1 and maintain a 
lookup table for the tail probability of for different values of the noise powers and path-loss 
exponents using (fTTj) . The lookup table is presented for a homogeneous 2-D SCS in Figure 0] as a 
plot of P ({t^v > l}) a g amst noise powers for different values of path-loss exponents. Further, 
in a homogeneous l-D SCS with a high BS density A , the equivalent SCS corresponds to a 
noise power N' — > according to Corollary [71 in which case, j, where —> corresponds 

to the convergence in distribution. Hence, in an interference-limited system (large Ao), the signal 
quality is measured in terms of j. In the following section, we study the effect of shadow fading 
on the j and at the MS in the SCS. 

V. Shadow fading in a SCS 

In Il25ll . (291, we have shown that the y at the MS is unaffected by shadow fading in a 
homogeneous 2-D SCS. Here, we present a theorem for obtaining the tail probability of j and 
given by (J2J), for any SCS. Again, without loss of generality, we consider the 1-D SCS. 

Theorem 3. When shadow fading in the form of i.i.d non-negative random factors, {^i}, is 
introduced to 1-D SCS with BS density function A(r), such that the ^i's are independent of the 
underlying Poisson point process of the 1 -D SCS, the resulting system is equivalent to another 



1-D SCS with a different BS density A(r) = ^ ?r 



Such an equivalence is valid as 



long as 



< oo. 
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Proof: The proof is provided in Appendix |Cj ■ 

Further, the j and at the MS in the 1-D SCS with shadow fading is the same as 

that in the equivalent 1-D SCS given in Theorem [3] where there is no shadow fading. Hence, 

using Theorem [TJ and Corollary [6] on this equivalent 1-D SCS, we can find $/ C n-i (lo) and 

vt) 

$/ c n.-i (uj), respectively, and by plugging these into (J7]) and (fT7|) , the tail probability of j 
and are computed. The following remarks are an immediate consequence of Theorem [3] 

Remark 3. When shadow fading in the form of i.i.d non-negative random factors, {^i}, is 
introduced to a homogeneous l-D SCS with BS density A , the resulting system is also a 



homogeneous l-D SCS with BS density AqE^, 



as long as 



< oo. 



Remark 4. In a homogeneous l-D SCS with BS density A , 

(a) shadow fading has no effect on the j at the MS (From Remark [3] and Remark Q] ), and, 



(b) the effect of shadow fading is completely captured in the noise power term of the 
(From Remark [3] and Corollary [7]). 

Figure [2] shows that the tail probability of j at the MS in a homogeneous l-D SCS with and 
without shadow fading are identical. Remark |4] provides a thoretical proof for this observation. 
Now, consider the following interesting example concerning log-normal shadow fading in a 
homogeneous l-D. 

Example 1. Consider a homogeneous 2-D SCS with an average BS density A , where each 
BS is affected by an i.i.d log-normal shadow fading factor which in dB has a mean and 
standard deviation a. Using Remark [31 the equivalent homogeneous 2-D SCS has a BS density 
A = A exp 0^r~\ Note that A > A , V a, e, and from Remark HI the introduction of shadow 
fading actually improves the signal quality, measured in terms of j^y. 

In the following section, we present a theorem that helps generalize the results for any arbitrary 
path-loss model. 

VI. Generalization to Arbitrary Path Loss Model 

So far, we have assumed the path loss model to be For some situations (like indoor 
channels), the path loss model could be different. However, the results in prior sections can still 
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be applied if we transform an arbitrary path loss model to the canonical path loss model ~ by 
the following theorem. 

Theorem 4. Let h(x), x > 0, be a monotonically increasing function with derivative h'(x) > 0, 
V x > and inverse function the 1-D SCS with BS density function A(r) and path loss 

model is equivalent to the 1-D SCS with a BS density function A(r) = ^ h -i^ cmd the 
canonical path loss model -| . 

Proof: Let R = h(R) be the equivalent BS location. Using the Mapping Theorem in ||27l , 
BS with locations R is also a Poisson point process, whose density is obtained below. For 
any non-homogeneous 1-D Poisson point process, E[N(t + s) — N (t)] = f* +s \{z)dz is the 
expected number of occurrences in the interval (t,t + s). Thus, 

E [N (r + s) - N (r)] = E [Number of BSs with Re (r,r + a)] 

= E [Number of BSs with R e (V 1 (r) , h' 1 (r + a))] 

= r ir+ "x {z)dz = r^'f'w'd, 

Jz=h-i(r) Jz=r h'{h 1 (Z)) 

Hence, the 1-D SCS with path-loss model and a BS density function A(r) is equivalent to 
the 1-D SCS with path-loss model — and BS density function A (r). ■ 

Once the new BS density is found, all results in prior sections apply after setting e — 1. For 
example, a homogeneous l-D SCS with path loss model R~ e is equivalent to the 1-D SCS with 
path loss model and BS density function, A'(r) = jr" 1 . This explains why $^ c ^-i (u) in 
Corollary [3] is only a function of ^ — 1. 

In the following section, we present several applications of the theory developed so far, in 
studying important scenarios in wireless communications pertaining to cellular system, adhoc 
systems, multiple access techniques, cognitive radio technology and overlay of wireless networks. 

VII. Applications in wireless communications 

Here, we study several scenarios in wireless communications by modeling the underlying 
system as a homogeneous l-D SCS with BS density A , where I = 1 models the BS deployment 
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along the highway, I = 2 models the planar deployments, and I = 3 corresponds to deployments 
in metropolitan and high-density residential areas. 

Cellular Communications: So far, we have studied the signal quality at the MS in a cellular 
system where the signals from the BSs are affected by arbitrary path-loss and shadow fading. 
Several other situations besides these can also be studied. Firstly, transmission powers and the 
antenna gains, {Ki}°Z v of all the BSs in the cellular system are seldom the same, in which case 
{K i ]^L l can be modelled as i.i.d. random variables with an appropriate distribution, independent 
of the underlying Poisson point process. By treating i^'s as shadow fading factors, their effect 
on the signal quality at the MS can be found. Moreover, all the remarks in Section [V] hold for 
Ki's as well. For example, consider the case where each BS has an ideal sectorized antenna 
with gain G and beam-width 9, such that BS's antenna faces the MS with probability j-, in 
which case K { = G, and otherwise Kj = 0. In this case, from Theorem [31 A = XqG^-^ is the 
BS density of the equivalent homogeneous l-D SCS. 

Multiple Access Techniques: Next, the signal quality at the MS in a cellular system employing 
different multiple access techniques can be studied. For example, in a code division multiple 
access (CDMA) system, the goal is to maintain a constant voice signal quality at the MS, 
which is done by power control. This goal is acheived by having the serving BS increase its 
transmission power by a = y (j) 1 , where a is the power control factor or the processing 
gain, y is the instantaneous signal quality at the MS, and y is the desired constant signal quality. 
In this formulation, a for each BS is a random variable and in general, the en's of nearby BSs 
are correlated. But if the correlation is small, the characteristic function of ( j) 1 computed in 
Section IIII-AI enables radio designers to approximately model the power needs to communicate 
with a MS in a SCS. In another formulation, if a is a constant factor by which the power 
of the serving BS is improved, its effect on the tail probability — at the MS is obtained by 
straightforward manipulations as P({«y > y \ e, /}) = K§_ (^) £ if y > a. 

Now, consider frequency division multiple access (FDMA) and time division multiple access 
(TDMA) based cellular systems. Let the available spectrum (in frequency for FDMA and in 
time-slots for TDMA) be divided into N parts called channel reuse groups (CG), and indexed as 
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k = 1,2, • • • ,7V, . Then, each BS is assigned one of the N CGs, such that the k th CG is assigned 
with the probability p^. In such a system, the MS chooses a CG that corresponds to the best ~; 
the BS in the CG that corresponds to the strongest received power is the desired BS, and the MS 
chooses it as the serving BS. The j at the MS in such a SCS is of interest to us. Note that this 
homogeneous l-D SCS is equivalent to N independent homogeneous l-D SCSs with constant BS 
densities Aq.Pi, • • • , XoPn, by the properties of Poisson point processes. The tail probability of 
j at the MS in such a system is given by P ({ £ > y\ e, N}) = 1 - [l - P ({ £ > y\ e})] N , 
where the tail probability on the right hand side is computed using Corollary [3] and (jTJ). 

Cognitive Radios: In cognitive radio technology, the cognitive radio devices (or secondary 
users) opportunistically operate in licensed frequency bands occupied by television (TV) transmitter- 
receiver pairs and wireless microphone systems (collectively called primary users). The inter- 
ference caused by secondary user transmissions is harmful for primary users operation, and is 
not acceptable beyond certain limits. Studying the nature of these interferences and formulating 
methods for addressing them has been an active area of research. The results in this paper are 
a rich source of mathematical tools for such studies. In Il30ll . we have extensively applied the 
results developed here, towards understanding the role of cooperation between the secondary 
users in ensuring that the interference caused by the secondary users are within the acceptable 
limits. The secondary users are modeled analogous to BS placement in homogeneous l-D and 
2-D SCS, and the tail probability of y at the primary user is characterized. Further, in the 
context of radio environment map (REM, ll30l and references therein]), we have highlighted the 
practical significance of the study of l-D SCS. A more detailed account of the study of cognitive 
radio networks incorporating the results from this paper is a topic for future research. 

Overlay Networks: The modern cellular communication network is a complex overlay of 
heterogeneous networks, such as macrocells, microcells, picocells and femtocells. This complex 
overlay network is seldom studied as is, due to the analytical intractabilities. In 0, OTTl . Il32ll . 
cellular systems consisting of macrocell and femtocell networks are analyzed. Using the results 
in our paper, the cumulative effect of all the networks constituting the overlay network, on the 
signal quality at the MS can be studied. A detailed study on this is set aside as a future work. 
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VIII. Conclusions 

In this paper, we comprehensively study the characterization of the j and at the MS in the 
SCS where the SCS is defined as a cellular system where the BS deployment in a given region 
is according to a non-homogeneous Poisson point process. Analytical expressions for the tail 
probabilities of - and j-^f at the MS are obtained for 1-D, 2-D and 3-D SCSs, where the 1-D, 2- 
D and 3-D SCS are mathematical models for BS deployments along the highway (1-D), in planar 
regions (2-D) and in metropolitan and residential areas (3-D), respectively. Further, a closed form 
expression for the tail probabilitiy of j is derived for the homogeneous cases of 1-D, 2-D and 3- 
D SCS. Several mathematical tools are developed, that help in characterizing the instantaneous 
y and j| ^ at the MS when the MS is moving in the SCS. These results are applicable for 
general shadow fading distrubutions and arbitrary path-loss models. This makes the results 
directly applicable for analyzing many different wireless applications which are characterized by 
uncoordinated deployments. The application of this theory has been demonstrated in the study 
of the impact of cooperation between cognitive radios in the low power primary user detection 
and can be found in 11301 . Future work will further explore the applications of the SCS model 
in the context of indoor femtocells, cognitive radios, and multi-tier or overlay networks. 

Appendix 

A. Useful Properties of 1-D SCS 

In this section, some of the frequently used properties of the non-homogeneous Poisson point 
process in the context of 1-D SCS with BS density function A(r) are listed. Note that all the 
proofs are not provided here, and can be found in 112211 . [|27ll . or in any of the elementary books 
on stochastic processes and Poisson processes. 

Fact 1. If Ri denotes the separation of the serving BS from the MS, then, the p.d.f. of R\, 

fRi (n) , is f Rl (n) = A (n) exp (- J^ Q A (s) ds) , V n > 0. 

Fact 2. If the BSs in the 1-D SCS are sorted in an ascending order of their distances from 
the MS denoted by random variables R 1 ,R 2 ,- • •, then, the distance of the (k — l) th nearest 
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BS from the MS is smaller in the usual stochastic order (Definition |4] in Section IIII-A1I) than 
the nearest k th BS, i.e. Ru-i <st Rk, cmd the p.d.f. of Rk conditioned on Rk-i is given by 
fR h \R*-i (rfc|r*_i) = A (r k ) exp(- f s t Tk i \{s)ds^ , V r k > r k -i, k > 1. 

Fact 3. When the separation of the k th closest BS, Rk = r k , the entire 1 -D plane may be divided 
into 2 regions, namely, [0,rk}and (r k ,oo). Any event in [0, r&] is independent of any event in 
(rjt, oo). Also, any event in non-overlapping intervals within the same region are independent of 
each other. The above statements follow from Definition [2] 

Lemma 3. Let {R k }T = i r ^P r ^sent the set of distances of BSs from the MS (indexed in the 
ascending order of the distance), Dk+\ = Rk+i — Rk be the distance between two adjacent BSs, 
and fD k+1 \R k (d\ r; A(s)) be the p.d.f. of D k+ i conditioned on R k = r corresponding to the 1-D 
SCS with BS density A(s), s > 0. Then, 

f Dk+1 \ Rk {d\r;X(s)) = e-^ +d ^\( r + d), and (19) 



faD k+1 \aR k (d' \r' ] A(s)) = f Dk+1 \R k [ d' r'] -A(^) J 

\ Hi Hi J 



(20) 



Proof: The p.d.f. in 4H follows from FactEl Further, f aDk+1 \ a R k (d'\r'; A(s)) = \f D \ Rk (£| r -] 

r'+d! 

= e sr ±A r r_±«Lj = e Jr> a \ a j ±^ , where (a) is obtained by expressing 

the p.d.f. in terms of D k+ i conditioned on R k , (6) follows from (jT9j) , (c) is obtained by suitably 
manipulating the integral in the exponential, and finally, fl20|) is obtained from (c). ■ 

B. Proof of Theorem [7] 

We first derive the expression for the characteristic function of Pj conditioned on Ri, <£> p J |^ 1 (u\ 
and use this to obtain the characteristic function of (~) \ Q^cy 1 (^)- 

Derivation for ^Pj^ (u\ri): Consider a region [0, Rb) in the 1-D plane such that Rb 3> r±, 
and parition into PI : [0, ri] and P2 : (ri,R B ). Next, partition F2 into iV non-overlapping 
intervals P2j : (r^ = r\ + (j — 1) Ar, r\,- + Ar = n + jAr) , where j = 1, 2, 3, • • • , N and 
Ar = 3azri. As N — > oo, Ar — > and it becomes unlikely for more than one BSs to exist 
in any partition P2j. If X, denotes the contribution to Pi by BSs in the j th partition, then Xj 
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is a Bernoulli fa) random variable defined as P ({Xj = Kr~ £ + o (Ar) | R\ = ri}) — Pj, and 
P ({Xj = 0| R\ = r±}) = 1 — py, where pj = A(r.,)Ar + o(Ar) is the probability that there is 
exactly one BS in P2j, j = 1, 2, • • • , N. The characteristic of Xj is 

$^1* (w|ri) = (A fa) Ar + o (Ar)) e K^T +0 (Ar)) + x _ A Ar + (Ar) 
fbl / ■ j?- — e \ f c i (iujKrJ 6 ) 1 



k + X fa) Ar (e iuKr i " - l) + o(Ar) ( = ) l + A fa) Ar ^ ^ -i-L + o(Ar) 



vfc=l 



(«0 



[iuKr- £ ) 



k ' 



exp ( A (r,) Ar ( £ ^ ) ) + o(Ar), (21) 



where (a) follows from the definition of the characteristic function, (b) is obtained by approxi- 
mating the sum of two little-o functions with an equivalent little-o function, (c) is obtained by a 
Taylor's series expansion on the exponential function, and (d) is obtained by adding the necessary 
terms to form an exponential function and compensating for them in the little-o function. 

Next, as Rb —> oo, the total interference power is Pj = J2f=iXj, where {Xj}™ =l is a set 
of independent random variables, independent of R\ due to our construction of PI, P2, and 
P2j,j = 1, 2, ■ • ■ and Fact [3] in Appendix lAl The characteristic function of Pj given R\ = r\ is 

N 

®Pi\Ri(u\ri) = lim lim TT^x i^ (w|ri) 

Rb — >-oo iV— >oo 

Ar-tO j=2 

N / / oo ( ■ ts -e\ k " 



^Too exp ( E ( A fa) Ar [ £ ) ) ) + °( A -) 



Ar->0 

-Kb 



i=2 V \fc=i 



R \im exp Qf B A (r) (V^ r ~ £ - l) dr^j 

lim exp lim / A (r) > ; dr 

R B ->oo \m->oo J r=ri ^ n! x 

/ A faKf [ R ° _ n£ , ' 

exp lim > ; — lim / Alrir dr 



-R E 



Notice that Qpj^ (ufa) exists if and only if lim J = B A(r)r ne dr < oo, V n > 1 and V r x > 0, 



in which case $ Pl \ Rl (ufa) = exp \ J2n=i IZ n X ^ r ™ £rfr 
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Derivation for (oj): Firstly, recall that j = where P s = KR X e , and $ , £r ii 



ri). Hence, $/ C n-i (cj) 



fying the expression for $ 



E 



Ri 



Ri 



and (J4j) is obatined by simpli- 



-Ri 



ri). Lastly, A(r) should satisfy JJ =0 A(s)rfs < oo, for the 



p.d.f. of i?i in Fact [2] of Appendix [A] to exist. 

Remark. The expression for the characteristic function can be derived from the characteristic 
functional of the Poisson point process in [27, Section 3.3, Page 31], but the proof technique 
presented here is emphasized, as it will be used in solving similar problems in the future. 



C. Proof of Theorem \3\ 

Let R = , where R is the random variable representing the distance from the MS to a BS 

in the 1-D SCS with a BS density function A (r), ^ is the shadow fading factor corresponding to 
the BS, and R is the corresponding equivalent distance. Using the product space representation 
and Marking Theorem in Il27ll . R also corresponds to the 1-D SCS with a BS density function 
derived below. 

For any 1-D SCS, E [N (t + s) - N (t)], the expected number of BSs in the interval (t, t + s) 
is called the mean function and is given by E [N (t + s) - N (t)] = J t +S X(z)dz. Let E[Number 
of BS in R E (r, r + s)] be denoted by E [N (r + s) - N (r)]. Then, 



E [N (r + s) - N (r)] = E 



(r+s)*J 



A (z) dz 



(b) 



(r+s) 



En 



where (a) is obtained by rewriting the expectation with respect to each realization of and (b) 



is obtained by exchanging the order of integration and expectation in (b) as Eq 
oo. Hence, R's corresponds to the 1-D SCS with a BS density function A (r) = Ey 



tye A {Z^ 



< 



D. Simulation Methods 

In this section, the details of simulating the SCS are presented. A single trial in simulating 
the BS placement for the 1-D SCS with BS density function A(r) in the region of interest which 
is a subset of the 1-D plane denoted by S, involves the following steps: 
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1) Generate a random number M ~ Poisson (J s A (s) (is), which is the number of BSs to be 
placed in S for the given trial. 

2) BS placement: For homogeneous 1-D SCS, generate M random numbers according to a 
uniform distribution in the range of S. If A(s) does not correspond to a homogeneous 1-D SCS, 
if A max = supA(s), then general a random number ?/ ~ Uniform(0, A max ) and another random 
number x according to a uniform distribution in the range of S. A BS is placed uniformly at x, 
only if y < Xq(x). This process is repeated until M BS are placed. 

3) Compute the received power at the MS for each BS using the path loss exponent e. The 
shadow fading in the SCS is incorporated by multiplying each of the received powers with i.i.d. 
random number generated according to the distribution of the shadow fading factor. Finally, j 
at the MS corresponding to this trial, is computed according to ©. 

For all the simulations in this paper T = 20000 trials are used unless specified otherwise. 
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Figure 1. Invariance of the performance of homogeneous l-D SCS on the BS density. 
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Figure 2. Comparison of Simulations with the analytical results 
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